Facilitated Asymmetric Exclusion 
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We introduce a class of facilitated asymmetric exclusion processes in which particles are pushed 
by neighbors from behind. For the simplest version in which a particle can hop to its vacant 
right neighbor only if its left neighbor is occupied, we determine the steady state current and the 
distribution of cluster sizes on a ring. We show that an initial density downstep develops into a 
rarefaction wave that can have a jump discontinuity at the leading edge, while an upstep results 
in a shock wave. This unexpected rarefaction wave discontinuity occurs generally for facilitated 
exclusion processes. 

PACS numbers: 02.50.-r, 05.40.-a 



In the asymmetric exclusion process (ASEP), sites of 
a lattice are occupied by single particles, each of which 
can hop at a fixed rate to a neighboring vacant site on 
the right [l|-l5| ■ This versatile model describes many sys- 
tems, including traffic [6[j9[, ionic conductors [lOj, and 
RNA transcription [ll|, [l^. Despite its simplicity, the 
properties of the ASEP are rich and deep. For exam- 
ple, a density that increases with x leads to a propagat- 
ing shock wave, similar to a traffic jam that propagates 
along a congested road. Conversely, when the initial den- 
sity drops quickly as a function of a;, a rarefaction wave 
arises in which the drop gradually smooths out, as occurs 
in stopped traffic after a stoplight turns green. Macro- 
scopic aspects of these phenomena can be understood 
from hydrodynamic theories [13J, [14 1, while the fluctu- 



ations about these macrostates continue to be actively 
investigated [l5|-tl9|. 



FIG. 1: Illustration of occupancy facilitated asymmetric ex- 
clusion. Particles that are eligible to hop to the right are dark, 
while immobile particles are shaded. This configuration con- 
tains islands of lengths 2, 3, and 1 (left to right). 

In this work, we investigate facilitated asymmetric ex- 
clusion. We primarily focus on occupancy facilitation in 
which a particle can hop to its vacant right neighbor only 
if its left neighbor is also occupied (Fig. [1]) . This model 
was proposed by Basu and Mohanty 20] in the context 
of non-equilibrium absorbing state phase transitions. We 
also investigate distance facilitation in which the rate at 
which a particle hops to a vacant right site is a decreasing 
function of the distance between a particle and its closest 
left neighbor. 

The notion of facilitated exclusion is part of a general 
class of ASEP models in which the hopping rate of a 
particle depends on more than just the occupancy of the 
neighboring site [6|-|8|, 120, l21| . For example, in glassy dy- 
namics the particle mobility decreases as the local density 
increases j22l | . Conversely, the presence of nearby parti- 



cles may increase hopping rates; for example, in molecu- 
lar motor models a moving particle can exert a hydrody- 
namic force that pushes other particles along [23| . More- 
over, a subset of phase space in occupancy facilitated 
exclusion can be mapped onto the ASEP of extended ob- 
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a model that was formulated to 



mimic the traffic of ribosomes along RNA. 

In occupancy facilitation, a mean-field hypothesis for 
the current is J — p^{l — p); the expression accounts for 
the presence of two particles and one vacancy and repre- 
sents a natural generalization of the current J — p{l — p) 
in the ASEP. As we show below, the current in facili- 
tated exclusion actually has a very different density de- 
pendence. We also develop a hydrodynamic description 
for an initial density step and predict that a rarefaction 
wave develops a discontinuity at the leading edge. Fi- 
nally, we provide a general criterion to understand this 
unexpected phenomenon in the framework of distance fa- 
cilitation. 

Finite Ring: We first determine the density dependence 
of the current on a finite ring in occupancy facilitation. 
The key to understanding the steady-state spatial distri- 
bution of particles is the notion of islands. An island is 
a string of occupied sites that are delimited at both ends 
by vacant sites (Fig. [1]). Each hopping event transforms 
a triplet • • o into • o •. Depending on the occupancy 
of the next site, the number of islands either increases, 
• •oo — > •o»o, or remains the same, ••o»— >»o»», 
but cannot decrease. Thus the system eventually reaches 
a state where the number of islands is maximal. 

For p < -k, the constraint that the number of islands 
can never decrease ensures that the system eventually 
reaches a static state that consists of immobile single- 
particle islands. The approach to the final state has a 
rich time dependence 27[, particularly in the marginal 
case p — k where the number of active particles asymp- 
totically decays as t~^/^ (see also Refs. 28|-|30[). 

In the p > I steady state, the requirement that the 
number of islands is maximal ensures that adjacent va- 
cancies must be separated by at least one particle. Fur- 



thermore, configurations that contain the maximal num- 
ber of islands are equiprobable. Indeed, let P{C) be the 
steady-state probability of being in a maximal- island con- 
figuration C. Then the stationarity condition is 

P{C)J2RiC^C')=J2P{C')RiC'^C), (1) 
c c 

where R{C —>■ C) is the evolution rate from configura- 
tion C to C" . Since _R = 1 if an evolution step is allowed 
and otherwise, we need to count the number of ways 
into and out of a configuration to solve Eq. ([T]). 

The evolution out of a configuration is triggered by 
triplets of the form • • o at the right edge of any island 
of length > 2. The system can evolve into the configura- 
tion C from another maximal-island configuration by the 
process ••o»— )-»o»». This evolution can only happen 
at the left edge of an island of length > 2. Hence there 
are an equal number of terms on both sides of Eq. ([T|) . 
If P{G) are equal for all configurations, Eq. ([T|) is clearly 
satisfied. Thus, all maximum island configurations are 
equiprobable in the steady state. 
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FIG. 2: Illustration of number of places that V vacancies can 
be placed among A'^ particles (filled circles) with: (a) site i 
occupied or (b) vacant. 

The probability of a maximum-island configuration 
therefore equals C^^, where C is the total number of such 
configurations with N particles and V vacancies on a ring 
oi L = N -\-V sites. To determine C, consider an arbi- 
trary site that we label by i. If this site is occupied, there 
are N possible locations between the N particles to put 
the V vacancies (Fig. [2]). If site i is unoccupied, there 
are A^ — 1 possible places to put the remaining V — 1 
vacancies. In both cases, we cannot put more than one 
vacancy between consecutive particles or else the number 
of islands would not be maximal. The number of such 
configurations is therefore given by 



C 



N 
V - 



(2) 



To obtain the steady state current, consider the flow 
across a link between arbitrary adjacent sites i and i-\-l. 
For a particle to move across this link, the consecutive 
sites I — 1 and i must be occupied while site i -I- 1 must 
be vacant. We now enumerate the number of maximum- 
island configurations that are consistent with the pres- 
ence of this triplet by noting that there are N — 2 places 
between the remaining particles to place the V — 1 re- 
maining vacancies so that no two vacancies are adja- 
cent. Thus the number of allowed configurations con- 
sistent with the presence of this triplet is {yZi)- The 



J = 



iv-l) . (1-P)(2p-1) 



C 



(3) 



with p = j; held constant in the limit N^L — > oo. 
(This result can be mapped into an equivalent ex- 
pression for the current in the ASEP of extended ob- 
- 26| ; we return to this correspondence 
The current is zero at p = i, since the system 



jects [11, 

below.) itie current is zero at p 

eventually reaches the static state of alternating particles 

and vacancies. The current is also zero at p = 1 where no 

evolution is possible. The maximal current arises when 

P* = ^, where J{p*) = J„,ax = 3 - 2\/2 « 0.1716. 
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FIG. 3: Current versus density for occupancy facilitation. 
The smooth curve is the prediction ((Sjl, while the circles are 
simulation data from 10^ realizations on a ring of 10^ sites. 

We can also determine /„, the density of islands of 
length n. Using the same enumeration that gave the 
number of allowed configurations, there are V — 2 remain- 
ing vacancies that can be distributed among the N —n~\ 
places between the rest of the particles so that there are 
no consecutive vacant sites. There are ( y "2^^) such con- 
figurations. Since each configuration has equal weight, 
the density of islands of length n is 



(N-n-l\ 
\ V-2 ) 
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where the latter equality applies for n ^ L: both Eqs. ([3]) 
and ([ll were also derived in Refs. 20|, l2]J by indepen- 
dent methods. The island length distribution decays as 
A", with A = (2p — l)/p, rather than X = p that occurs 
for a random particle distribution. Since (2p— l)/p < p, 
long islands are suppressed compared to a random dis- 
tribution; this feature is a consequence of the constraint 
that the number of islands is maximal. From these is- 
land probabilities we recover the particle density from 
p = ^ nin , while the current J can alternatively be ex- 
pressed as the probability to have an island that contains 
at least two particles, J = X]n>2 ^n- 



Density Step: Let us now study the evolution of a den- 
sity step on the infinite fine by occupancy facilitation. 
Initially, the density to the left of the origin is p_ , while 
the density to the right is p+. For a downstep, where 
P- > p+, the density profile within a hydrodynamic de- 
scription evolves by the continuity equation -gf + §j =0, 
which we may solve by the method of characteristics [13| . 
The solution is a function of a scaled variable, z = x/t, so 
p{x,t) = f{z). Using the steady state current expression 
([3]) for the flux, we flnd that the scaled proflle is com- 
posed of distinct segments in which the density is either 
constant or given by / = (2 + z)^^/^. Thus the density 
profile is 



Different density profiles arise depending on whether 




z < z_ 
Z- < z < z+ 
z > z+ . 



(5) 



The position of the left interface 2_ is determined from 
continuity: (2 -|- z_)~^/^ = p_. When p_ > p* = -i=, 
we have z_ < 0. In this situation, the density at the 
origin p(0) is universal and it coincides with the density 
p* that maximizes the current in Eq. ([3])- Therefore the 
number of particles that penetrates into the region a; > 

is N{t) = J[p{Q)]t = J„,axt. 
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p+ < 2 or p+ > 



In the former case the right interface 

2 



is located at z-|_ = [2 — 3/9+ — 2-^/(1 — /0+)(l — 2p+)]/p-^. 
As z passes through z+ the density jumps from the value 
(2 + z+)~^/^ to p+. For example, when {p-,p+) = (1,0), 
z+ — J and the magnitude of the density drop is | 
(Fig. U). The discontinuity at the front of a rarefac- 
tion wave arises because the leading particle cannot move 
unless "pushed" by neighboring particles from behind. 
Consequently, the density at the leading edge must be 
non-zero. 

For p+ > -^j this jump discontinuity disappears, and 
the density profile is everywhere continuous. Continuity 
a.t z — z+ now gives p+ = (2-|-z_(_)~^", which manifestly 
solves Eq. (|6]). For this class of rarefaction waves, the 
density is sufficiently large ahead of the wave that the 
leading edge can get pulled ahead and there is no need for 
a pileup of particles from behind to push the wavefront 
forward. 

To study shock waves, we suppose that ^ < p- < p+ 
and consider a large region that includes the interface. 
The particle influx to this region is J_ , while the outflux 
is ,/+ . The net flux must equal the change in mass c{p- — 
p+) inside this region, where c is the shock wave speed. 
Hence c — (J_ — J+)/{p- — p+). Using the expression 
([3]) for the current, the shock wave speed is 



(P-P4 



(7) 



The shock propagates to the right if p_ < (2/3+)"^ and 
to the left otherwise. 

Our results can be extended to a more stringent occu- 
pancy facilitation in which r consecutive sites to the left 
of a particle must be occupied for a particle to hop to a 
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vacant right neighbor [2; 



(For example, for 



x/t 



the rightmost particle in o • • • o cannot move, while the 
update ••••o— >-»»»o»is possible.) A steady state 
with a non-vanishing current and a maximal number of 
islands, each of length > r, arises when p > j^. All 
such conflgurations are again equiprobable. 



FIG. 4: Scaled density profile of facilitated exclusion starting 
from the step initial condition p_ = 1 and p+ = 0. The simu- 
lation data is based on 10^ realizations for three representative 
times and is visibly indistinguishable from the prediction of 
Eq. dSl) wheni= 1.5^^ 



To locate the right interface z^, we apply the con- 
straint that the initial mass within [z_,z+] must equal 
the mass in this region at some later time plus the net 
influx into this region. In scaled units, this conservation 
statement is 



P-Z- +P+Z+ 



dz 



V2" 



+ J- 



■h 



(6) 



with J± — J{p±) 



time 
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FIG. 5: Equivalence between occupancy facilitated exclusion 
for the case r = 1 and the ASEP of dimers (rectangles) that 
hop to the left. 

We now discuss the connection between occupancy fa- 



cilitation and the ASEP of extended objects [11|, [1^ [21 



2J-|26|. In the maximal-island steady-state regime with 
density p > ^, we may equivalently view a particle fol- 
lowed by a vacancy as an extended object of length k — 2 
which hops to the left. Since vacancies cannot be adjacent 
in the steady state, these extended objects obey exclusion 



p- 




z < z_ 


I(z) 


Z— 


< z < z+ 


p+ 




z > z+ 



and perform a simple ASEP (Fig. [5]). This connection 
continues to hold with the more stringent r-tuple occu- 
pancy facilitation. In the region of phase space where 
p > ;^q:j- , the steady state behavior of the system maps 
to the ASEP of extended objects with length k = I + r. 
Finally, we treat distance facilitation. To find the cur- 
rent even for the simple example in which the hopping 
rate equals i~^ , where £ is the distance to the nearest 
left particle, is challenging. The hydrodynamic behavior, 
however, is robust and the rarefaction wave discontinuity 
always arises [27|]. We can demonstrate the universality 
of this phenomenon from basic features of the current- 
density relation. We know that J(0) = J(l) = and 
we expect that J{p) has a single maximum at some den- 
sity p*. Applying the scaling ansatz for the continuity 
equation shows that either p is constant or ^ = z. The 
rarefaction wave therefore has the form 



P(^) 



where I{z) is the inverse function of z — ^. Differ- 
entiating this relation with respect to z in the region 
Z- < z < z^ gives JppPz = 1. If Jpp is everywhere neg- 
ative (as in the standard ASEP), then pz must also be 
negative. Thus the density p{z) continuously decreases 
until it reaches p+. However, if Jpp is positive at some 
low density, then pz would become positive. Thus the 
smallest possible density Pmin in a rarefaction wave oc- 
curs at the point where Jpp vanishes. If Pnun > P+, there 
must be a jump discontinuity at the leading edge. 

Thus an inflection point in the current-density rela- 
tion signals a rarefaction wave discontinuity. Such an in- 
flection point must exist for any facilitation mechanism, 
since the Jpp{p*) < at the maximum p* and Jpp > 
for small p. One such example is J ~ p"+^ as p — >■ that 
arises for distance facilitation with hopping rate i~°'. 

In summary, facilitated asymmetric exclusion has fea- 
tures that are dramatically different from simple asym- 
metric exclusion. The most prominent is the jump dis- 
continuity at the leading edge of rarefaction waves. This 
phenomenon arises in a broad class of cooperative trans- 
port models with facilitated dynamics. 
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